function result = test_gp(train=10, sigma=0.01, theta=0.1, l=10)

load chirps.mat

X = chirps(1:train, 1);
t = chirps(1:train, 2);

mse = 0;
for i = 1:5
    x_star = chirps(10 + i, 1);
    result = gp(X, t, x_star, sigma, theta, l);
    real = chirps(10 + i, 2);
    F = result.F;
    mse += (real - F)*(real-F);
end
result = mse;
